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We obtain exact asymptotic results for the disorder averaged persistence of a Brownian particle 
moving in a biased Sinai landscape. We employ a new method that maps the problem of com- 
^1^' puting the persistence to the problem of finding the energy spectrum of a single particle quantum 

Hamiltonian, which can be subsequently found. Our method allows us analytical access to arbitrary 
values of the drift (bias), thus going beyond the previous methods which provide results only in the 
limit of vanishing drift. We show that on varying the drift the persistence displays a variety of rich 
asymptotic behaviors including, in particular, interesting qualitative changes at some special values 
_^ of the drift. 

O ■ PACS numbers: 05.40.-a, 46.65. +g, 02.50.-r 
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I. INTRODUCTION 

, Persistence, i.e., the probability that a fluctuating field does not change sign upto time t has been widely studied in 
recent years in the context of nonequilibrium systems [1]. A wide variety of results, both theoretical and experimental, 
are available for pure systems. In contrast, there have been very few studies of persistence in disordered systems. A 
notable exception is the study of persistence, theoretical [2] as well as numerical [3], in disordered Ising models. In 
'"^ , this paper we study analytically the persistence in another disordered system namely the celebrated Sinai model [4], 
^ ■ but in the presence of an additional arbitrary drift. We show that as one varies the drift parameter, the disorder 
O I averaged persistence displays a wide variety of rich behaviors which undergo qualitative changes at certain special 
values of the drift. 

The Sinai model [4] is perhaps one of the simplest models of disordered systems where various disorder averaged 
physical quantities exhibit rich and nontrivial behaviors and yet, can be computed analytically [5]. Thus the Sinai 
model serves the role of 'Ising' model in disordered systems. In this model a Brownian particle undergoes diffusion in 
J — . ■ presence of a random time-independent potential. The position x{t) of the particle evolves via the Langevin equation, 

O . 

rT\ ■ dx dU 

(N 

, where ri{t) is the thermal noise with {ri{t)) = and {■q{t)r]{t')) = 5{t — t') and U{x) is the external potential. In 
' the biased Sinai model one considers the potential to be simply U{x) = —^x + ^/aB(x) where B{x) represents the 
^ , trajectory of a Brownian motion in space, i.e., B{x) = £^{x')dx' with {S,{x)) = and (^(a;)^(x')) = b(x — x'). The 
y ' parameter /i represents the bias or the drift and a measures the strength of the disorder. Thus the particle is subjected, 
in addition to the thermal noise ri{t), an external position dependent random force F{x) — —dU/dx — /i+^(a;). Various 
(— ( • physical quantities in the Sinai model have been studied before [5]. Here our aim is to compute the persistence P{xq, t) 
Q , defined as the probability that the particle does not cross the origin upto time t starting at the initial position xq > 
O ' at i = 0. Evidently this quantity will depend on the realization of the underlying disorder potential and will, in 
general, vary from one sample of disorder to another. Our final goal is to compute the disorder averaged persistence 
P(xo, t) as a function of both xq and t for different values of the drift /i. 
^ . , Physicists have recently studied the persistence in the Sinai model using various methods which include an exact 
' probabilistic approach suited for unbiased Sinai model [6], a study of an equivalent lattice model with random 
hopping rates [7] and also by employing a real space renormalization group method [8] . All of these methods provide 
asymptotically exact results, but only in the limit of vanishing drift, i.e., when fi 0. Unfortunately, extension of 
these existing physical methods to extract explicit asymptotic results for arbitrary nonzero /i seems rather difficult. 
Mathematicians, on the other hand, have studied some aspects of a related quantity namely the distribution of the 
first-passage time in the Sinai model for nonzero drift and some rigorous results seem to exist [9-12]. However, these 
mathematical methods are extremely technical and difficult to follow. What is lacking, so far, is a unified physical 
approach which, besides reproducing the known results in a simple and transparent way, provides exact asymptotic 
results for all /i and yet simple and powerful enough to be easily generalizable to other problems. The purpose of 
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this paper is to provide such an approach. The heart of our approach hcs in mapping the problem of computing the 
persistence in the Sinai model with arbitrary drift to finding the spectrum of a single particle quantum Hamiltonian, 
which can subsequently be done exactly. 

Apart from presenting an unified approach valid for arbitrary drift, there arc two other physical motivations for 
this work. First, it is well known that the Sinai model displays a range of intereresting anomalous diffusion properties 
as one tunes the drift n through certain finite 'critical' values [5] . It is therefore theoretically interesting to know how 
the persistence behaviour changes as the drift is varied through these 'critical' values. Secondly and perhaps more 
importantly, the Sinai model with a nonzero drift has numerous physical applications [5] including the diffusion of 
electrons in disordered medium in presence of an electric field, glassy dynamics of dislocations in solids, dynamics 
of random field magnets, dynamics near the helix-coil transitions in heteroploymers. The most recent application 
of the biased Sinai model has been to understand the dynamics of denaturation of a single DNA molecule under an 
external force [13]. Persistence seems to be a natural quantity to study in these systems and hence we expect that 
the analytical results presented in this paper will be useful in many of the physical situations mentioned above. 

The paper is organized as follows. In section II, we present our general approach. A detailed discussion of the 
pure case {a = 0) with nonzero drift is presented in subsection II-A which will help us anticipate the general features 
of persistence in the disordered case studied later in subsection II-B where we illustrate the mapping to a quantum 
mechanics problem. In section III, we discuss the results for the disorder averaged persistence for positive drift (/U > 0). 
The results for the negative drift (/x < 0), fundamentally different from the positive drift case, are detailed in section 
IV. We conclude in Section V with a summary and outlook. The details of the derivation of the eigenvalue spectrum 
of the quantum Hamiltonian are presented in the appendix- A. In appendix-B, we present an alternative derivation of 
the disordered averaged persistence in the case of positive drift. The details of the second order perturbation theory 
for negative drift are presented in appendix-C. 

II. GENERAL APPROACH 

Consider the particle whose position x{t) evolves via the Langevin equation (1) starting initially at x{t = 0) = xq- 
The persistence P{xQ,t) is the probability that the particle does not cross the origin upto time t starting at xq. It 
is also useful to define the distribution of the first-passage time F{xo,t*) which is simply the probability that the 
particle hits the origin for the first time at t = starting initially at xq. The distribution F(a;o,t*) is related to 
the persistence P{xo,t) via the simple relation, P{xo,t) = 1 — Jq F(xo,t*)dt*. This follows from the fact that the 
integral Jq F(a;o,t*)rft, sums up the probabilities of all the events when the particle hits the origin before time t and 
when subtracted from 1 (which is the total probability), the resulting quantity, by definition, is the persistence. Our 
objective is to first compute F{xo,t^) or its Laplace transform and then use the above relationship to compute the 
persistence P{xo,t). 

In order to calculate the first-passage time distribution F{xo,t^.) we employ a powerful backward Fokker-Planck 
approach which has been used before to study the persistence in the unbiased Sinai model [6] as well as in other 

contexts [14]. It is instructive to start with a more general quantity, Qp{xo) = (e ^ ^[^(* where () denotes 

the thermal average, y[a;(t)] is an arbitrary functional and t» denotes the first-passage time, i.e., the time at which 
the particle first hits the origin starting initially at xq at t = 0. Note that if we choose V{x) = 1, then Qp{xo) = 
(e~P**)^^^ = e~P^* F{xq, t^)dt^ is simply the Laplace transform of the first-passage time distribution F{xo, i*). For 
convenience of notations, we will henceforth denote the initial position xq = x and the first-passage time t» = t. 

A differential equation for Qp{x) can be derived by evolving the particle from its initial position x over an infinites- 
imal time dt. This gives Qp{x) = ((1 — pVdt)Qp{x + dx)) where dx is the displacement of the particle in time dt 
from its initial position x. Using Eq. (1) one gets dx = F{x)dt -\- r]{0)dt where F{x) = —dU/dx = IJ. + ^/a^{x) is the 
random force. Expanding Qp{x + F{x)dt + ri{0)dt) to order dt and averaging over ri{0), one arrives at the backward 
Fokker-Planck equation, 

Since here we are interested only in the first-passage time distribution, we will henceforth set V{x) = 1. Note, however, 
that this method is powerful enough to deal with the statistical properties of any arbitrary functional V{x) of the 
stochastic process. We will also assume, without any loss of generality, that the initial position x > 0. For a; < 0, 
one will obtain the same results by changing the sign of the drift jj.. The equation (2) is supplemented with the two 
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boundary conditions: (i) Qp{x = 0) = 1, since if the particle starts at a; = 0, obviously its first-passage time t = 
and (ii) Qp{x oo) = 0, since the first-passage time t ^ oo if the particle starts at x = oo. 

In the next two subsections we discuss the solutions of the differential equation (2) respectively for the pure case 
{a = 0) and the disordered case (cr > 0). 



A. Pure case with nonzero drift 



It is instructive to discuss first the pure case with nonzero drift (fi ^ 0) in the absence of the random potential 

((7 = 0). The results for the pure case will help us anticipate what to expect for the disordered case later. Solving 

Eq. (2) with F{x) = /U we get Qp{x) = e"^''"^'' ^^^J"^ that satisfies the required boundary conditions. Since 
Qp{x) = e~P*F{x,t)dt, we need to invert the Laplace transform which gives 

F{x, t) = a;(27rt3)-V2e-(^+^*)V2t. (3) 
Using P{x,t) = 1 — F[x, t')dt', one then gets 

P{x, i) = 1 - ^ /* t'-'/'e-("+^*')'/2t'dt'. (4) 

Let us analyze what happens for large t in the three separate cases (i) > 0, (ii) /u < and (iii) = 0. 

(i) For positive bias away from the origin (/i > 0), the particle eventually escapes to oo with a nonzero probability 
and P{x,t) P{x) as t ^ oo. Taking t — » oo limit in Eq. (4) one easily obtains this eventual persistence 'profile' 

P{x) = 1 - 

(ii) In the opposite case (/i = — < 0), it follows from Eq. (4) that as t ^ oo, P{x,t) « 
y27^x|/x|-2r3/2e-(^-|''l*)'/2t. In this case, a more useful information is contained in the asymptotic first-passage 
distribution F{x,t). Prom the exact expression of F{x,t) in Eq. (3), we find that in the appropriate scal- 
ing limit X — > oo, i ^ oo but keeping x/t fixed, the first-passage distribution approaches a delta function, 
F{x,t) 6{t — x/\fx\). Equivalently the Laplace transform, Qp{x) which also follows directly from 



the expression of Qp{x) = exp — + 2p— x in the correct scaling limit x 0, p but keeping the 

product px fixed. Thus in this case, at late times, the particle essentially moves ballistically with velocity and 
crosses the origin for the first time at t = x/\iJ,\. 

(iii) In the unbiased case (/i = 0), we recover from Eq. (4) the well known exact result [15], P{x,t) = erf(a;/\/2t) 
where erf (a;) = ^ e~" du is the error function. 

In the next subsection, we switch on the quenched disorder (cr > 0) and examine the consequences on the asymptotic 
properties of the disorder averaged persistence . 



B. Disordered Case with nonzero drift 



Unlike the pure case, we can no longer solve the differential equation (2) exactly for the disordered case since 
F{x) = jj + \/o'£,{x) now has an x-dependent random part. To make further progress we first make a Hopf-Cole 
transformation, Qp{x) = exp [— Zp{x')dx'] . Clearly Zp{x) = —dlogQp{x)/dx is a slope variable. By construction, 
Qp{x) automatically satisfies the boimdary condition Qp{x = 0) = 1. Substituting this form of Qp{x) in Eq. (2), we 
find that the slope variable Zp{x) satisfies a first order stochastic Riccati equation 

= zlix') - 2 [m + V^a^')] zp{x') - 2p. (5) 

The right hand side of the above equation contains a multiplicative noise term and we will interpret it according to 
the Stratonovich prescription. Note that since Eq. (5) is a first order equation, Zp{x') at an arbitrary x' will be fully 
determined as a functional of the noise history {i,{x')}, at least in principle, provided the value of Zp is known at 
some 'initial' point. Note that this 'initial' point can be anywhere. The program would then be to substitute this 



3 



fully determined functional to evaluate the integral Jq Zp{x')dx' and thereby determine Qp{x) = exp [— Zp(.x')(i.x'] . 

Subsequently one would perform the disorder average Qp{x) where the overbar indicates an average over the noise 
history (x')} for fixed x. 

However, there is one problem in implementing this program namely the 'initial' value of Zp{x') is not specified. 
Consequently the solution of the first order equation (5) will involve an unknown parameter, i.e., the 'initial' value 
of Zp{x'). There exists, however, a rather nice trick to get around this diSiculty. This trick has been used before 
in the Sinai model in various contexts [5,6,20]. It is useful to outline this trick in the present context. To use this 
trick, we first fix x in the definition Qp{x) = exTp[—w{x)] where we have defined w{x) = Zp{x')dx' . Keeping x 
fixed we then make a change of variable, t = x — x'. Thus, when a;' — > oo, r — > — oo and when a;' — > 0, r — > a;. 
Besides, x' = x corresponds to r = (see Fig. 1). Then w{x) = Zp{x')dx' 
Zp{t) = Zp{x — r). In this new r variable, the equation (5) becomes, 

dZp{T) 



■ Zp{x — T)dT = Zp{T)dT where 



dT 



zI{t)+2 fl+^^{T) Zp{T)+2p. 



(6) 



where ^(r) = ^{x — r) and r G [— oC',.t]. Note that (^(r)) = and {^{t)^{t') = 6{t — t'). To simplify further, we 
substitute Zp{T) = exp [(^(r)] in Eq. (6) and find that the variable (^(r) satisfies a much simplified stochastic equation 
containing only additive noise (and no multiplicative noise) 



^ = b{cj>) + 2V^^~(t), 



where the source term h{(j)) is given by 



(7) 



(8) 



What did we gain in the change of variable from x' G [0, oo] to r e [— oo,a;]? The point is that we have, via this 
change of variable, gotten around the problem of 'initialization' of the original variable Zp{x') at any 'initial' point. 
To see this, let us consider the Eq. (7) which is valid in the regime r G [— oo,a;]. We can interpret this equation 
now as a simple Langevin equation describing the evolution of the position (P{t) of a classical particle with 'time' r 
starting from r = — oo. Of course, we still do not know the value of (/)(t) at t = — oo. The point, however, is that 
this initial condition at r = — oo is completely irrelevant. No matter what this initial condition at r = — oo is, it is 
clear from Eq. (7) that eventually when r is far away from its starting point r = — oo, the system will approach a 
stationary state. This is because the Eq. (7) describes the noisy 'thermal' motion of a particle in a classical potential 
Uc\{<l>) = — Jq b{u)du = e'^ — 2iJ,(j)+2pe~'^ — {2p+ 1). Hence the particle will eventually reach the equilibrium and the 
stationary probablity distribution of (j) is simply given by the Gibbs measure. 









' I r<t> 


Pst (<?!)) =^ exp 




= A exp 


- b{u)du 



(9) 



where 6(0) is given by Eq. (8) and A is a normalization constant such that Pst{(f>)d(j) = 1. For later purposes we 
also define Pst(</') = 4'o{4') where 



V'o (<?!>) = V^exp 



1 

4^ 



b(u)du 



(10) 



the function b{(j)) is given by Eq. (8) and A is such that ^Q{<j))d<p = 1. 

So now we know that starting at r = — oo with arbitrary initial condition, by the time the system reaches r = 0, 
it has already achieved the stationary measure. But our task is not yet complete. We now have to evolve the system 
via its equation of motion (7) from t = to r = .t (knowing that at r = the distribution of 4> is given by the Gibbs 
measure in Eq. (9)) and evaluate the disorder average 



Qp{x) = E [exp[-w(a;)]] = E 



exp[ 



Jo 



dr] 



(11) 



where w{t) = Jq e'^^'^ ^dr' as defined earlier and E[..] denotes the expectation value of the random variable inside 
the parenthesis. Let us introduce the quantity i?(0, r) = E^[e^^'^^'^^] which denotes the expectation value e"'''"' at 
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time T with (/"(t) = (/). More precisely, if Pj[w{t) = w,(f){T) = (f>,T] denotes the joint probabiUty distribution of the 
variables u'(t) and 0(t) at time r, then R{(j), t) = J e~'^'^Pj[w, (p, rjdw. We have introduced the additional parameter 
A for later convenience whose value will be eventually set to A = 1. Note that if we set A = 0, R{(l), r) is simply the 
probability distribution of (j) at time r. Thus from now on, we will refer to the A = case as the 'free' problem. When 
A = 1, it is clear that, 



Qp{x) = E [exp[-«;(a;)]] 



R{(j>-, x)d(j}. 



(12) 



The advantage for this small detour in introducing the new quantity R{4>, r) is that one can now write down a 



Fokker-Planck equation for R{(t), r 
R{cI),t) = [exp[-A/;e^(^')rfr']' 
equation, 



in a straightforward manner. In fact, incrementing r to r + dr in the definition 
and using the Langevin equation (7), we find that r) satisfies the following 



dR d^R 



dR 



[h'{(fy+\e^] R, 



(13) 



where b'{(p) = db/d(p with b{(f)) given from Eq. (8). Note that at t = 0, w{Q) = and hence R{cf),0) is just the 
probability distribution of which is given by the Gibbs measure R{(j),Q) = Pst(0) in Eq. (9). Starting with this 
initial condition at r = 0, we need to evolve the equ ation ( 13) upto t = x, determine R{(l>, x) and then integrate over 
(f) in Eq. (12) to finally obtain the desired quantity Qp(x). Note that for A = 0, the Eq. (13) is simply the ordinary 
Fokker-Planck equation for the probability distribution of (f> in the 'free' problem. 

We next substitute R{4>,t) = exp ^ b{u)du G{(p,T) in Eq. (13) to get rid of the first derivative term on the 

right hand side of Eq. (13) and find the following evolution equation for the Green's function G{^, r). 



dG 

a7 



2a 



|62(</,) + l6'(</.) + Ae^ 



G, 



(14) 



where a = 1/4(7 and G{(j),0) = exp —-^ b{u)du R{(j),0) = \/A'ipo{(j)) , using the Gibbs measure in Eq. (9). To 
solve Eq. (14) we make the standard eigenvalue decomposition 

G(</.,T)=^CB5i^(<^)e-^"''^ (15) 



where the eigenfunctions gsifp) satisfy the Schrodinger equation 

1 d^9E{^) 



+ 



'^b\<P) + ^b'{<P)+aXet' 



gE{4>) = EgEicj}), 



(16) 



with 6(0) given by Eq. (8) and a = l/4c7. The coefficients c^'s in Eq. (15) are determined from the initial condition, 
G{(p,0) = VAipoicp)- Using orthogonality of eigenfunctions one finds. 



/oo 
d(j)g*E{4>)^l)Q{(j)) = y/A< gsli^o 
-OO 



>, 



(17) 



where we have used the standard bra-ket notation of quantum mechanics. Note also that if we consider the 'free' 

problem by sotting A = 0, it is easy to verify from Eq. (16) that there is an cigcnfunction with energy E = which 
corresponds to the stationary state of the 'free' problem. This zero energy eigenfunction is given precisely by tpo{<P) 
in Eq. (10) and the Gibbs measure is just the square of this eigenfunction, Pst{4') = i^oW- 

Substituting the c^'s from Eq. (17) in the decomposition equation (15) and setting finally r = a; we get the following 
expression of R{<p, x) in terms of the eigenfunctions. 



i?(0, x) = ^ < gE\i)o > gE{4')e 



4a Ex 



(18) 



By integrating R{(j),x) in Eq. (18) over cp, we finally obtain the disorder average Qp{x) in a compact form. 
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R{(l), x)d(f) = I d(j)'4}l{(t))gE{(})) 5^ < 9E\i)o > e" 

-cso J — oo ^ 



ia-Ex 



-iaEx 



= <^o|e-^"^"|Vo>, 
where the quantum Hamiltonian H in the </> basis is given by, 



H 



(19) 



(20) 



Here Hq 



1 
2 



+ ^b'{(t>) is the Hamiltonian of the 'free' problem (corresponding to A = 0), Hi = Xae'^ 



is like a perturbation Hamiltonian and b((j)) 



2ii + 2pe' 



The exact formula Qp{x) = dte~P^F{x, t) =< ipo\e~'^'^^^\'tpo > in Eq. (19) is, in fact, the central result of this 
paper. This result tells us that the Laplace transform of the disorder averaged first-passage time can, in principle, 
be fully computed for arbitrary starting position x, arbitrary p (and hence for arbitrary t) and also for any value of 
the drift provided one can compute all the eigenvalues and the corresponding eigenfunctions of the Hamiltonian H 
in Eq. (20). In other words, the calculation of the disorder averaged persistence is reduced to finding the spectrum 
of the quantum Hamiltonian H. In the next two sections we show how this spectrum can be determined in limiting 
cases which lead to exact asymptotic results (large t limit) for the persistence P{x, t) for any arbitrary drift /x. 



III. EXPLICIT RESULTS FOR POSTIVE DRIFT (/x > 0) 

In this section we focus on the positive drift (/x > 0) case. We have seen in Section II-A that for the pure case, 
the particle eventually escapes to infinity with a nonzero probability when there is a positive drift (/x > 0) away from 
the origin. This escape probability P{x) is precisely the persistence P{x,t) in the limit t oo. Due to a nontrivial 
dependence of this probability P{x) on the initial position x, we call P{x) the persistence profile. Note that from 
the relationship, P{x,t) = 1 — F{x,t^,)dt^ where F{x,t) is the first-passage time distribution, it follows that the 
persistence profile is given by P{x) = 1 — Jg°° F{x, t^)dt^ = 1 — Qo{x) where we recall that Qp{x) = Jq°° e~^*F{x, t)dt 
is just the Laplace transform of the first-passage time distribution. 

In the disordered case with positive drift, one would expect a similar behavior namely for each sample of disorder, the 
particle will eventually escape to infinity with a nonzero sample dependent probability P{x). The disorder averaged 
persistence profile is then given by P{x) = 1~ Qq{x). Using this relationship one can then compute, in the large t 
limit, the exact time-independent persistence profile by setting p = in the general formula for Qp{x) in Eq. (19) 
derived in Section II-B. For p = 0, we get from Eq. (8), &(</>) = — e"^ -I- 2/x. Substituting this b{(j)) in Eq. (20), setting 
A = 1 and simplifying, we find H = Hm + where v = [ija and Hm is a generalized Morse Hamiltonian given by, 

2d<^^ 2^ 2 ^ ■ ^^^> 

It then follows from Eq. (19) that 

= e-""""^/^ < Vo|e-^"*""|Vo > . (22) 

To evaluate the matrix element in Eq. (22) explicitly we need to know the spectrum of the generalized Morse 
Hamiltonian Hm- Fortunately this spectrum can be fully determined. This calculation is done in details in appendix- 
A. Here we just summarize this spectrum and use the results to compute P{x) = 1 — Qo{x) explicitly. 

The spectrum of Hm consists of two parts: a discrete part with negative energies that correspond to the bound 
states and a continuous part with positive energies corresponding to the scattering states (see appendix- A). The 
nature of the spectrum depends on the parameter v = fi/a. It turns out that there is a critical value I'c = 2 such 
that for 1/ > Pc, the spectrum has both the bound states and the scattering states. In contrast, for i/ < Vc, there 
are no bound states and only scattering states exist. We notice that a similar behaviour was obtained in the study 
of transport properties of the Sinai model [20]. The eigenvalues and the corresponding eigenfunctions are given as 
follows. 
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Bound States: The bound states are labelled by an interger n. The eigenvalues are given by 

Er, = -\\vl2 - 1 - nf, n = 0, 1, . . . ^2 - 1] 



(23) 



where [m] indicates the integer part of m. Clearly this discrete spectrum exists provided v > 2. The corresponding 
normalized eigenfunctions are given by 



9n{4>) = bne-^^^W^ (2ae*) 



(24) 



where Wa,/3{x) is the Whittaker function [17]. The normalization constant 6„ can also be computed exactly (see 
appendix A) 



2 _ 2(7(i^ - 2 - 2n)T{n + 1) 



r{u-l- n) 



(25) 



where r(x) is the standard Gamma function. 



Scattering States: The scattering states have positive energies labelled by the wavevector q, Eg = g^/2 with 
< g < 00. The corresponding eigenfunctions are given by 



= b{q)e-'^/'W^.^ (2ae*) , 



where the coefficient b{q) is given by 



1 „, r{l-./2-iq) 
Ti-2iq) ■ 



(26) 



(27) 



This coefficient b{q) is chosen such that in the limit (j) ^ -co (where the quantum potential in the Hamiltonian Hm 
in Eq. (21) vanishes) the eigenfunction gq{4>) approaches a plane wave form, i.e., gq{4>) — > [e*^''* + r(g)e~'''^] as 

(f) — oc', where e"^"^ represents an incident wave travelling in the direction of positive (j) and e"*'''^ represents the 
reflected wave travelling in the opposite direction with r{q) being the reflection coefficient (for details see appendix- A). 

Having obtained the full spectrum of Hm we are now ready to compute the persistence profile P{x) = 1 — Qo{x). 
Expanding the right hand side of Eq. (22) in the energy basis of Hm and using the results on the spectrum of Hm 
summarized above, we get 



Qo{x) 



-crv'^x/2 



[v/2-l\ 



/ -I roo 

V I < 5n|^o > |2e2-(-/2-i-")^- + / dq\< 5,|Vo > ?e 

n=0 



2 -2q'ax 



(28) 



a result which is valid for all x and for all > 0. The function Vo('^) is already known. In fact, for p = we find 
from Eqs. (8) and (10) the following normalized expression 



M4>) 



exp 



4(7 2^ 



(29) 



Using this expression of V'o(9!') and the eigenfunctions in Eqs. (24) and (26) one can easily evaluate the matrix elements 
< gn\i^o > and < fiiglV'o >• For the bound states we get 



/oo 
-oo 



V2^r>) 

where 6„ is given by Eq. (25). Similarly for the scattering states we obtain 



T{n + l)T{u -l-n), 



< 9q\^0 >- 



d(l)ipo{(j))gg{(j)) 



b{q) 



=r {v/2 - iq) r {v/2 + iq) , 



y2rfX^ 

with b{q) given by Eq. (27). Substituting these matrix elements in Eq. (28) we get our final expression. 



(30) 



(31) 
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Qoix) 



[v/2-l] 



J2 &^r2(n+l)r2(:/- l-n)e2'^("/2-i-«)=x^ / dq\b{q f \r {i^/2 - iq) fe'^"''''' 

n=0 Jo 



(32) 



where 6„ and bg are given respectively by Eqs. (25) and (27). Substituting Eq. (32) in the relation P{x) = 1 — Qoix) 
then gives us the exact persistence profile valid for any a; > and any /i > 0. 

It is instructive to derive explicitly the tails of this profile P{x) for small x and large x. Consider first the limit 
a; — > 0. While we can use the general solution in Eq. (32) to derive the small x behavior, it is easier to consider the 
original equation (19) which for small x gives Qo{x) ^ 1 — 4:ax < ipolHlipQ >. Since H = Hq + Hi and moreover 
since tpo{(l)) is a zero energy eigenfunction of Ho, we get Qq{x) ^ 1 — iax < ipolHilfpo >■ Expanding the matrix 
element in the ^-basis, using Hi = ae'^ (setting A = 1 in Eq. (20)) and the expression of ipoift^) from Eq. (29) and 
evaluating the resulting integral, we get Qq{x) ^ 1 — 2ah'x. Using ly = /i/cr, we get P{x) = 2 fix as x ^ 0. Thus we 
obtain an interesting result that the slope 2/i characterizing the linear growth of the profile near a; = is completely 
independent of the disorder strength a. In fact, the small x behavior of the disorder averaged persistence profile is 
identical to that of the pure case. 

We now turn to the other limit x ^ oo. Here we use Eq. (32). First consider the case when i/ > 2. Then we know 
from Eq. (32) that there exist bound states. In that case it is evident that for large x, the term corresponding to the 
lowest energy boimd state (n = 0) will be the most dominant term on the right hand side of Eq. (32) . Retaining only 
this leading n = term in Eq. (32) and using 6q = 2a/r(i^ — 2) we get, Qo{x) = E^e~'^'^'^^~^'>^ as a; — * oo for > 2. 
Consider now the opposite case when v < 2. In this case there are no bound states and there is no contribution 
from the discrete sum on the right hand side of Eq. (32). The only contribution is from the integral representing the 
scattering states. For large x, the most dominant constribution to the integral will come from the small q regime. 
Expanding the Gamma functions for small q, we find after preliminary algebra, Qo{x) = Ai,{2ax)~^^'^e~'' '^^/^ as 
X ^ oo for < 2 where Ai, is a constant (see below). Exactly at = 2, we get from Eq. (32), Qo{x) = e~'^'^^ /\/2nax 
for large x. Let us summarize the three different types of large x behaviors of the persistence profile, 

> 2, 



1 - P{x) i 



7^.'-'^^ - = 2' (33) 



<"^/2 V <2, 



where A^, = 7r'^/^r^(i//2)/r(i/) [1 — cos(7ri/)]. Evidently the shape of the profile in Eq. (33) for large x changes as v 
varies through the critical point Vc = 2. The reason for the existence of this critical point is evident from our analysis. 
Essentially it happens due to the loss of bound states as v decreases from > 2 to < 2. Note that this critical 
behavior at finite u = Vc = 2 could not be derived by the RSRG method. The RSRG method is valid only in the 
1/ — > limit where the exact result in Eq. (33) coincides with the RSRG results [8] . 

We conclude this section by pointing out that it is possible to have an alternative derivation of the disorder averaged 
persistence profile P{x) for /x > by a completely different method. This method relies on mapping the calculation 
of the persistence profile to calculating the disorder average of the ratio of two partition functions in the Sinai model. 
This mapping makes use of certain mathematical properties of the Brownian motion. The average of this ratio of 
partition functions was already computed before [9,18] in a different context. Using those results one can then recover 
the results in Eq. (33). The derivation of this mapping is presented in appendix-B. 



IV. EXPLICIT RESULTS FOR NEGATIVE DRIFT (/x < 0) 

We now turn to the the case when /x < which corresponds to a drift towards the origin since the initial position 
a; > 0. The situation here is very different from the positive drift /i > case discussed in the previous section. For 

< 0, we expect from the analogy to the pure case that for each sample, the particle will definitely cross the origin 
as f ^ oo no matter what the starting position x is. Hence iov ji < the persistence P{x,t) ^ as t ^ oo for 
all X, unlike the /U > case where the persistence approaches a time independent profile P{x,t) — * P{x) as t ^ oo. 
Therefore what is interesting in the /i < case is to compute the asymptotic behavior of P{x, t) for large t. Recalling 
the definitions P{x,t) = 1 — F{x,t»)dt» and Qp{x) = e~P*F{x,t)dt where F{x,t) is the first-passage time 
distribution, we see that the analysis of the large but finite t limit of the disorder averaged persistence P{x, t) requires 
an analysis of the p ^ limit of the Laplace transform Qp{x) [rather than exactly at p = where Qo{x) = 1 trivially 
for /X < 0]. In fact we will see later in this section that in this limit, P{x,t) or equivalently F{x,t) display a variety 
of scaling behaviors as one tunes the relevant parameter u' = —i^/a. 
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As in the case of /i > 0, the starting point of our analysis for /i < is the central result in Eq. (19) which is valid 
for all iJ.. Unlike the /j. > case, we can not, however, put p = straightway in the Hamiltonian H in Eq. (20). 
To derive the time dependent asymptotics for large t, we need to analyze the spectrum of H keeping p small but 
nonzero. Unfortunately, for nonzero p, it is hard to determine the full spectrum of H exactly. Fortunately, however, 
it is possible to extract the leading asymptotic behavior as a; ^ oo and p ^ without too much trouble. To see this, 
we consider the energy eigenvalue decomposition in the second line on the right hand side of Eq. (19) and find that 
as a; — > oo, the leading contribution comes from the ground state energy £^0 of H, 

Q^^\<H9o>fe-'^^°^ (34) 

where \go > is the ground state of H. It is clear from Eq. (34) that to evaluate the large x asymptotics we just need to 
compute the ground state energy Eq of the Hamiltonian H in Eq. (20). Unfortunately, even the ground state energy 
Eq is hard to compute exactly for arbitrary p. One can, however, make progress in the p — > limit. To see this we recall 
that \tpo > in Eq. (34) is the exact zero energy eigenstate of the 'free' part of the Hamiltonian Hq (see the discussion 
after Eq. (17)), i.e., Ho\ipo >— 0. Knowing this exact fact, we can then determine the ground state energy Eq of the 
full Hamiltonian H = Hq + Hi by treating Hi = Xae'^ (where eventually we will set A = 1) as a perturbation to the 
'free' Hamiltonain Hq. For example, to first order in the perturbation term, we get Eq = 0+ < i/'ol^^i IV'o > where 
indicates the fact the ground state energy of the unperturbed Hamiltionain Hq is exactly 0. We then decompose this 
matrix element < Vo|-ffi|V'o > in the (f) basis and use Hi = ae'^ (setting A = 1) to obtain Eq = a 'tpQ{4>)e'^d4>. The 
normalized wavefunction tpo{^) is already known from Eq. (10) and has the following explicit expression. 



/ i^p)"'/^ 



e'^ i^' p ^ 
"4^ ~ 2"*^" 2^^ 



(35) 



where Ka{x) is the modified Bessel function of index a [17] and we have used the definition f' = —n/a. Using this 
explicit form of il)o{(j>) and carrying out the integral in Eq = a ipQ{<j))e'^d<j) we get, to first order in perturbation 
theory, 

^Ki^{V2p/a) 

The result in Eq. (36) is only upto first order in Hi. It is not clear, so far, why one should stop at first order only. 
In other words, we have not yet specified in what sense Hi is 'small' compared to Hq. Note that in deriving Eq. (36) 
we have not yet taken the p limit. In the limit p 0, using the asymptotic properties of Bessel functions [17] in 
Eq. (36), we find 

Eq ^ \ -\p\ogp v' = 1, (37) 
[ aB^.p"' < z/' < 1, 

where we recall that a = l/4cr and the constant B^/ is given by 

2^-'''T{1- ly') 

B^' = ^2.'-ir(^.) • (38) 

Thus, in general, Eq ^ p" as p — > where a = 1 for v' > I and a = for i/' < I with additional logarithmic 
corrections at t^' = 1. Hence basically Eq is 'small' for small p. Naively one would expect that if, indeed, we carry 
out the perturbation theory in Hi to higher orders, the resulting terms will be lower order in p as p ^ 0. This naive 
expectation, fortunately, turns out to be true. In fact we show in appendix C how to estimate the second order term 
for small p and it turns out to be at least of 0(p"+^) and hence negligible compared to the first order term (~ 0(p")) 
in the p ^ limit. For small p, this argument therefore justifies in keeping only the first order term in the perturbation 
theory in evaluating Eq. Note that the eigenfunction \go > also gets modified from the 'free' eigengunction [^'o > due 
to the perturbation term. 

Substituting the p ^ results from Eq. (37) in Eq. (34) for large x and using a = l/4cr, we then get three different 
types of scaling behaviors depending on v', 

-pxla{v' - 1)] v' > 1, 

plog(p)x/(j] i^' = l, (39) 
-B^^p^'x] v' <l, 
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where B^/ is given by Eq. (38). Note that to leading order in smaU p, we need to keep only the zcroth order term in 
the amplitude | < ^o\9o > P of the exponential in Eq. (34). To the zeroth order l^o >= IV'o > and hence to leading 
order this amplitude is exactly 1 since \ipo > is normalized. Thus it is evident from Eq. (39) that for v' > 1, the 
correct scaling limit is x ^ oo, p ^ but keeping the product px fixed. On the other hand, for v' < 1, the correct 
scaling combination is p'^ x. 

It turns out that the exact asymptotic results in Eq. (39) were also derived recently by mathematicians by using 

completely different methods which involved rather heavy mathematical machineries. The first derivation is due to 
Kawazu and Tanaka who used the so called Kotani's formula and Krein's theory of strings [10,11]. However their 
method didn't permit the explicit evaluation of the constant B^/. More recently, Hu et. al. [12] presented yet another 
completely different derivation by mapping the persistence problem with negative drift onto that of a Bessel process 
and then using some theorems on this Bessel process. Hu et. al. managed to compute the coefficient B,^' explicitly. 
Note, however, that the constant B^/ in ref. [12] has an apparently rather different looking form than our expression 
in Eq. (38) and it requires a bit or work to show that indeed they are exactly identical. 

To derive the asymptotic properties of the first-passage time distribution F{x, t) in the time domain, we need to 
invert the Laplace transform in Eq. (39) with respect to p. For v' > 1, the Laplace inversion is trivial and we get 
F{x, t) = 5 [t — xj a(v' — 1)) in the scaling limit x oo, t — > cx3 but keeping the ratio xjt fixed. This result is very 
similar to the pure case. The delta function indicates that at late times the particle basically moves ballistically with 
an effective velocity a{y' — 1). Similarly for v' = 1, by inverting the Laplace tranform we find that in the scaling limit 
X — > oo, t oo but keeping the ratio xlog{x)/t fixed, F(x,t) = S{t — x\og(x)/a). The situation becomes somewhat 
different for i/' < 1. In this case the Laplace inversion indicates that in the scaling limit a; ^ oo, t ^ oo but keeping 
the ratio x^^'^ /t fixed, the first-passage time distribution approaches a scaling form F{x, t) ~ ^f{t/x^^" ). The scaling 
function f{y) can be formally written in terms of Levy function L,yi{y) and we get f{y) = yL^i{y). The Levy function 
is formally defined by the Bromwich integral (see the appendix of ref. [5]), 



1 
2^1 



d-\-ioo 



d—i 



(40) 



d line in the complex s plane lies to the right of all the singularities (branch cuts) but is 



where d is such that the s 
otherwise arbitrary. 

An explicit expression of the Levy function can bo obtained only for a few special values of v' [5] 
1/2, we get f{y) = yL,/2{y) = f-^/^y/ " 



For example. 



for i^' = 1/2, we get f{y) = 2/^1/2(2/) = e ^/^^/^/27rj/. This indicates that F{x,t) — xe ^ '^*/v27rt'^ and hence the 
persistence P{x,t) = 1 — F{x,t')dt' = eif{x/V2i). We thus arrive at an amazing result that for v' = 1/2, i.e., 
H = —(t/2, the disorder averaged persistence has the same asymptotic behavior as the pure case without drift [ji = 0, 
(7 = 0) (as derived in section II- A)! A similar coincidence at this special value of i'' = 1/2 was also noted in the 
context of occupation time distribution in the Sinai model [20]. Thus fi = —a/2 seems to be a special line in the 
{fi — a) plane where the Sinai model with drift shares the same asymptotic properties as the pure unbiased Brownian 
motion. Another solvable point is v' = 1/3, where we get f{y) = yLi^^{y) = ^"^ Ki/^ (o^o/y/y) with the constant 



r(2/3) 
r(i/3) 



3/2 



For general v' < I, while we can not calculate the scaling function f{y) explicitly, the behaviors at the tails can 
be easily determined. For example, first consider the limit y = tjx^l^ 00, i.e., when t » x^l^ . Using the 
large y behavior of the Levy function [5], Li,'{y) w Bi,'r{l + u')sin{'Ku')/'Ky'' we get F{x,t) « pox/t" +^ where 
/3o = u'2^-''' /T{u')a'^'''-'^. The persistence P{x,t) = 1 - F{x,t')dt' then behaves in this limit as 



F{x,t) 



■2v' 



t» X 



(41) 



^ for large t where the persistence exponent 9 = u'. We now turn to the other 
tjx^l^ 0, i.e., when t << x^^^ . Using the properties of the Levy function 
g-(i-..')/^'C/^27r(l-j/')C where C = [y"' /u'B^, 



indicating a power law decay P{x, t) t 

tail of the scaling function f(y) when y = 

near y ^0 [5], we get f{y) = yL^-{y) ~ e-^^-"' ^l"' ^ / s/2tt{1 - v')C where ( = [y"' /u' B^-Y'^'--'' K This indicates an 
essential singularity at y = 0. Using this asymptotic behavior of F{x, t) in the relation P{x, t) = 1 — J* F{x, t')dt' we 
find the following behavior for the persistence. 



P{x,t) ^ 1 - , ^ (^) 



-1/2(1-!.') 



exp 



l-u' 



1/(1-,.') 



t «x 



(42) 
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where /3i = {ly' B,,,)^/'-^-"'^ and B^, is given by Eq. (38). 

Let us summarize the main behavior of the disorder averaged persistence for i/' < 1. We find that for i/' < 1 , there 
exists a single time scale tx ^ x^^" depending on the initial position x. For t » t^, the persistence P{x,t) decays 
as a power law with an exponent 9 = v' and the amplitude of the power law depends on x as in Eq. (41). In the 
opposite limit when t « tx, the persistence drops extremely sharply from its initial value P{x,Q) = 1 as indicated 
by the essential singularity in the second term on the right hand side of Eq. (42) when t « x^l^ . 

Let us conclude this section on negative drift by one final remark. We note that even though we have assumed 
throughout this section that is strictly negative, we can safely take the limit /x ^ 0~ in Eq. (36) which gives 
£^0 ~ — l/[21ogp] in the limit ^) — > 0. It then follows from Eq. (34) that in the limit of vani shing d rift, one gets the 
asymtotic result Qp(x) ^ 1 + -^—x + . . . when x << — 1/logp. From this it follows that P{x,t) « 2cra;/Iogf for 
logt >> X, thus recovering the standard Sinai model behavior in the zero drift limit [8,6]. 



V. SUMMARY AND PERSPECTIVES 



In summary, we have obtained exact asymptotic results for the disorder averaged persistence in the Sinai model 
with an arbitrary drift. Our method maps exactly the problem of computing the persistence to the problem of finding 
the eigenvalue spectrum of a single particle quantum Hamiltonian. We have shown that it is possible to find this 
spectrum in certain cases which allowed us to obtain exact asymptotic results for arbitrary drift. We note that these 
results could not have been obtained from the existing physical methods (e.g. the RSRG method) which provide 
exact results only in the limit of zero drift. Our results show that there is a rich variety of asymptotic behaviors in 
the persistence as one tunes the drift. In particular, the asymptotics undergo interesting 'phase transitions' at certain 
critical values of the control parameter (the relative strength of the drift over the disorder), e.g., at = 2, = 
and u = —1. It would be interesting to extend the exact method presented here to calculate other properties in the 
Sinai model with finite drift, such as the persistence of a thermally averaged trajectory for which the results in the 
zero drift limit are already known [8,7]. 

We thank D.S. Dean for useful discussions. 



APPENDIX A: DERIVATION OF THE SPECTRUM OF THE HAMILTONIAN 



In this appendix we derive the spectrum of the generalized Morse Hamiltonian Hm given by Eq. (21). We show that 
the spectrum has a discrete part with negative energies which correspond to the bound states and also a continuous 
part with positive energies corresponding to scattering states. The eigenfunctions satisfy the Schrodinger equation. 



2 d(i2 



gE{(t>) = EgE{(t>)- 



(Al) 



Let lis remind the readers that a = l/4cr. It turns out to be convenient to make a change of variable y = 2aeJ 



Furthermore, let us substitute gsi^') 
equation. 



-0/2 



f{2ae''') in Eq. (Al). Then the function f{y) satisfies the differential 



+ 



1 1/4 + 2E' 



m = 0, 



(A2) 



where we have suppressed the E dependence of the function f{y) for notational convenience. We next consider the 
negative and the positive part of the spectrum separately. 



Bound States: The bound states are in the negative energy part of the spectrum. Let us substitute E = —7^/2 in 

Eq. (A2) where 7 is the eigenvalue to be determined. The second order differential equation (A2) with E = —7^/2 
is known to have two linearly independent solutions W,^-i ^{y) and W_i^ j{~y) where Wa,j3{x) is the Whittaker 
function [17]. Thus the most general solution of Eq. (A2) can be written as 

f{y) = D,W^^^{y) + D2W_._^^^{-y), (A3) 

where Di and D2 are arbitrary constants. For large argument, the Whittaker function is known to have the asymp- 
totic behavior [17], Wa,ii{x) ~ x"-e~^l'^. On the other hand the bound states must be normalizable and hence the 
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cigcnfunction gsi't') must vanish as <f> ^ ±00. The vanishing boundary condition at </) 00 indicates that the constant 
D2 = 0. Note that here we have assumed 7 > 0. If 7 < 0, then this boundary condition instead sets Di = 0. However, 
the resulting solution is the same. In other words the eigenfunctions corresponding to 7 and —7 are the same and not 
linearly independent of each other. Thus without any loss of generality we can assume 7 > and set D2 = 0. Thus 
the eigenfunction in terms of the original variable is given by 

gE{<l>) = D,e-^/^W^^^{2ae'^). (A4) 

Note that the eigenvalue 7 is yet to be determined. This is done by employing the vanishing boundary condition at 
the other tail, namely gE{4>) ^ a.s (p ^ —00. Using the small x behavior of the Whittaker function, Wa,i3{x) 
_^2^__^-/3+i/2 as x ^ in Eq. (A4) as ^ -00 we get 

Thus the eigenfunction diverges exponentially as (j) ^ —00. The only way the eigenfunction can satisfy the boundary 
condition gE{—oo) = is if the denominator r(l — i//2 — 7) in Eq. (A5) is infinite. This happens when the argument 
of the Gamma fimction is a negative integer 1 — v /2 — 7 = —n with n = 0, 1, . . .. This thus fixes the eigenvalue 
J = v/2 — 1 — n with n = 0, 1, . . .. However note that the condition 7 > indicates that the maximal allowed value 
for n is [i'/2 — 1] where [x] denotes the integer part of x. Thus finally the bound states have discrete eigenvalues 
E = —7^/2 = —{y/2 — 1 — n)^/2 and the corresponding eigenfunctions, labelled by n, are given by 

9n{<l>) = bne-'^/^W._^ ._^_^ {2ae'^) , (A6) 

where 6„ — Di is to be fixed from the normalization condition, J^^gn{4')d(p — 1- To perform this integral, we first 
use the fact that for positive integer n, one can rewrite the Whittaker function in terms of the Laguerre polynomials 
[17] and then use the following identity (see the appendix of [19]) 



One obtains 



2 2a{u-2-2n)T{n + l) 
= r(.-l-n) • ^^^^ 



Scattering States: We now turn to the positive energy part of the specrum and set E = g^/2 in Eq. (A2). 
The resulting differential equation, once again, has two linearly independent solutions Wu-i .^(y) and W_t^ iq{~y)- 

Note that the function W <^-i i^i^y) = W i^-i ^^{—2aeJ') ^ exp(e'^) as </) 00. Since the eigenfunction gE{<l>), even 
though non-normalizable for scattering states, can not diverge superexponentially a.s (f) —^ 00, this second solution is 
not allowed. Keeping only the first solution we get the eigenfuntions, now labelled by q, 

gM = b{q)e-^''W._^^,^ (2ae^) . (A9) 

The question is how to determine the constant h{q) in Eq. (A9). This is because, unlike the bound states, the 
eigenfunctions in Eq. (A9) are non-normalizable. To see this let us examine the behavior of near the tail 

— > — 00, as in the discrete case. Using the asymptotic properties of the Whittaker function we find that as ^ — > —00, 



gq{(t>) Ki) 



T{l-v/2-iq)^^''' ^ ^ T{l-v/2 + iq)^^''' 



(AlO) 



Clearly the functions gq{4>yii arc non-normalizable. Moreover, unlike in the case of bound states where the boundary 
condition at ^ = — 00 decides the discrete eigenvalues, in this case we have no such condition indicating that all 
possible values of 5 > are allowed. Note that, as in the discrete case, g > and q < correspond to the same 
eigenfunction and hence the allowed values of q lie in the range < g < 00. To determine the constant b{q) we 
note that in the tail (j) —00, the quantum potential in Eq. (Al) vanishes. The resulting differential equation with 
E = q^/2 allows plane wave solutions of the form 
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5,(.^)«-=[e'«^ + r(g)e-^«^], (All) 



where e'^'^ represents the incoming wave coming from </> = —oo and e ^'^'^ represents the reflected wave going towards 
(j) = —00 with r{q) being the reflection coeSicient. The ampHtude l/\/2n ensures that the plane waves ipq{x) — ■^==e*'^^ 
are properly ortho-normalized in the sense that < ipql^pq >= S{q—q') where 5{z) is the Dirac delta function. Comparing 
Eqs. (AlO) and (All) in the regime cj) — oo, we determine the constant b{q) up to a phase as 

This completes the derivation of the spectrum of the Hamiltonian Hm- 



APPENDIX B: ALTERNATIVE DERIVATION OF THE PERSISTENCE PROFILE FOR /x > 

For each sample of the disorder, the persistence proflle P{x) is related to the Laplace transform Qp{x) of the first- 
passage time distribution at p = via the relation, P{x) = 1 — Qq{x). The quantity Qo{x) can be obtained exactly 
by solving the Eq. (2) with p = 0, 

Q,{x) = l-p^ (Bl) 

where Zfj^{x) = e^u{x ^-^c partition function in a finite box of size x with U{x) = — F{x')dx' = ~^x + 

^JaB{x) being the random potential. It turns out to be useful to rewrite Eq. (Bl) in a slightly different form using a 
well known property of the Brownian motion: If B{x) is a Brownian motion, then B{x) — B{x') = B(x — x') where 
B{x) is another independent Brownian motion and = indicates that the random variables on both sides have the 
identical distribution, though they are not equal. Using this property and after a few steps of elementary algebra, 
one can rewrite Eq. (Bl) as, 

go(x) = 1 - P{x) = f ^ ' (B2) 

where .Z_^(a;) = e~'^^^'^'^^^^">dx' . Interestingly, exactly the same ratio as in Eq. (B2) has appeared earlier in 
other contexts and its average (over disorder) is known exactly [9,18]. Using these known results and setting v = /i/a, 
we get exactly the same asymptotic (large x) persistence profile as in Eq. (33) which was obtained in Section III by 
a completely different method. 



APPENDIX C: ESTIMATION OF THE ENERGY CHANGE DUE TO THE SECOND ORDER TERM IN 

PERTURBATION THEORY FOR < 

In this appendix we consider the case < and provide an estimate of second order contribution A£^2 to the 

ground state energy E = + AE'i + AE2 of the full Hamiltonian H — Ho + Hi given by Eq. (20), treating Hi = ae"^ 
(setting A = 1) as a perturbation to the unperturbed Hamiltonian Hq = —\-§^ + VQi'P)- The quantum potential 
Vq(0) is given explicitly by 



ae'^'l' + iv' - 1)6"^ + 4a f - p] - 2p(u' + 1)6-^ + Aap'^e'^''' 
\16a^ / 



(CI) 



In writing the explicit form of the potential we have substituted the expression of &(</>) from Eq. (8) in Eq. (20) and 
used the definition a = l/4cr. Let us recall that the ground state energy of the unperturbed Hamiltonian Hq is 
and the ground state wavefunction is '0o('/') given explicitly by Eq. (35) is section- IV. The first order contribution 
AEi =< tpo\Hit{jo > was already evaluated exactly in Eq. (36) and was to shown to scale as ~ as p ^ with 
a = 1 for z^' > 1 and a = u' for v' < 1. The goal of this appendix is to show that the second order contribution AE2 
is neglible compared to AEi as p — > 0. 
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Note that for p > 0, it is clear from Eq. (CI) that Vq((^) oo as (j) ^ ±00. This indicates that for p > the 
spectrum of Hq is discrete and consists of bound states only. Let ipn's denote the discrete energy eigenfunctions of Hq 
with corresponding eigenvalues denoted by e„. The second order contribution to the ground state AE2 then follows 
from the standard quantum mechanics, 



AEn 



n>0 



< V'o|gi|V'„ > 

eo - e„ 



= -E 



< Vol^fllV-n > 



(C2) 



n>0 



where we have used the fact that eo = as discussed earlier. Using the fact ei < 62 < 63 ... in Eq. (C2), one can 
immediately obtain an upper bound to —AE2, 



-AE2 < — ^1 < i^o\Hl\i^n > 
^1 n>0 



(C3) 



By adding and subtracting the n = term to the sum on the right hand side of the above inequality and using the 
completeness of eigenfunctions we get 



-AE2 < — 

ei 



< ^|JQ\Hl\'^l)o > -| < V'o|i?i|V'o > 



(C4) 



The quantity 5 =< ipo\Hi\tJjQ > — | < ^/jq | -ff 1 1 V-'o > P inside the parenthesis on the right hand side of the in- 
equality (C4) can be evaluated exactly. In general, for any m, one can express the matrix element < V'ol-f^i^lV'o >= 
a"* i)Q{(l))e'^'^ d(j} in the (j) basis by using Hi = ae'^. We then substitute the explicit expression of ipo{^) from Eq. 
(35). The resulting integral can be performed exactly using the identity [17] 



Jo 



7 



1//2 



(C5) 



In fact, throughout this paper, we have heavily used this identity. We then get < ipo\Hi^\ipo >= 
a'^{2p)"^/'^Kjn-u' /K^t (\/2p/cr). The expression of S requires the results for m = 2 and m = 1 and we 

get 



S{p) = la^v 



(C6) 



Using the expansion of Bessel functions for small arguments [17], it is easy to show that in the limit p — > 0, S{p) ^ 
for v' > 2, S{p) -p^ \ogp for 1/ = 2 and S{p) ~ p^' for 1/ < 2. 

Having established the behavior of S for small p, we now need to estimate the gap ei (the energy of the first excited 
state of Hq) for small p on the right hand side of the inequality —AE2 < S{p)/ei in Eq. (C4). To estimate the gap, 
we examine the quantum potential in Eq. (CI). It is convenient first to make a change of variable z = e'^ so that 
< ^ < 00. In this new variable the quantum potential in Eq. (CI) reads. 



Vq{z) = 



az 



+ - l)z + 4a 



16a2 



-p - 



2p{y' 



1) 4ap2 



(C7) 



The shape of this potential is shown in Fig. 2. Note that the potential Vq{z) has a minimum at 2 = Zo, determined 
from the equation dVQ{z)/dz = which gives, 2azQ + {v' —1)zq+2p{v' +l)zo—d,ap'^ = 0. In the limit p ^ 0, the only real 
root of this equation is at Zq k, Aap/{v' + 1). The value of the potential at this minimum, Vq{zq) —{2i/' + l)/8 + 0{p) 
as p — > 0. We also need to estimate the typical width of the potential W{e) at an energy e for small p (see Fig. 2). 
The points z±{e) where Vq{z) = e can be easily estimated for small p since in this limit one just needs to solve a 
quadratic equation and we get. 



z±{e) 



4a 



[u' + lT V^u' + 1 + 2e] 



p + 0(p2 



Hence the typical width of the potential at energy e scales as W{e) « z+{e) — Z-{e) 
finite level e. For small p, one can approximate the potential Vq{z) around its 



(C8) 

p in the p ^ limit at any 
^; = zo by a harmonic 
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oscillator potential, Vq{z) w — (2i/' + l)/8 + uj^z^ j2 where the frequency uj is estimated from the typical width , i.e., 
(^2^2^2 ^ 0(1). Since W(e) ~ -p, we find that the frequency scales as w ~ 1/p as p ^ 0. One knows that the gap 
between the first excited state and the ground state in a harmonic potential scales as e\ w. Thus we estimate that 
the energy of the first excited state scales as ei ^ ^ 1/p in the limit p ^ 0. 

Substituting this estimate of the gap in the inequality, —l\Ei < S{p)/ei and using the small p estimates of S{p) 
derived earlier, we find that as p — > 0, —AE2 < p^ for u' > 2, —AE2 < —p^logp for i/ = 2 and —AE2 < p" """^ for 
i/' < 1. Comparing these results with the first order contribution where AEi ^ p for v' > I and AEi ^ p'^ for v' < 1, 
we conclude that the second order contribution is negligible compared to the first order term in the limit p —^ 0. 
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X 







FIG. 1. The change of variable t = x — x' tor fixed x. The new variable t increases from r = — oo (when x' = oo) to r ; 
(when x' = 0) through the point t = (when x' = x). 
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-0.5 
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FIG. 2. The shape of the potential Vq{z) in Eq. (C7) is shown for parameter values: a = 1, v' = 0.5 and p = 0.1. The 
potential has a minimum at « = ~ P and a typical width W ~ pin the limit p — » 0. 
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